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Abstract 

We provide the first black hole solutions with Lifshitz asymptotics found in string theory. 
These are expected to be dual to models enjoying anisotropic scale invariance with dynamical 
exponent z = 2 at finite temperature. We employ a consistent truncation of type IIB 
supergravity to four dimensions with an arbitrary 5-dimensional Einstein manifold times a 
circle as internal geometry. New interesting features are found that significantly differ from 
previous results in phenomenological models. In particular, small black holes are shown to be 
thermodynamically unstable, analogously to the usual AdS-Schwarzschild black holes, and 
extremality is never reached. This signals a possible Hawking-Page like phase transition at 
low temperatures. 



1 Introduction 



The targets of the gauge/gravity duality have been vastly increased with the appreciation 
that it can be applied to strongly coupled condensed matter systems Pfl El E] ■ Among the 
phenomena of interest we encounter phase transitions governed by fixed points enjoying the 
anisotropic scaling t — > X z t while x — > A x. The parameter z is termed dynamical exponent, 
and z = 1 corresponds to the familiar scale invariance of the conformal group. Clearly, a 
departure from this value breaks conformality and thus the gravity dual cannot be an AdS 
space. 

Two different holographic realizations of this symmetry have been proposed. The first one 
is expected to model cold atoms at unitarity and the corresponding metrics were constructed 
in jU [5] . For the particular value z = 2 the isometry group of this space coincides with the 
invariance group of the Schrodinger equation and in consequence they are called Schrodinger 
geometries. 

The second avatar, Lifshitz geometries, is the focus of this paper. This class of metrics 
are named after the Lifshitz field theory that is supposed to arise as the long wavelength 
(continuous) limit of lattice models of strongly correlated electrons. They were introduced 
in [B] precisely as gravity duals of such systems and read 



where L sets the scale of the spacetime and for z = 1 one recovers AdS in Poincare coordi- 
nates. In addition to anisotropic scaling if r — > A -1 r, the isometries of these metrics include 
spacetime translations along with space rotations. 

Often the approach to condensed matter holography is phenomenological, with no man- 
ifest connection to string or M-theory. The matter content of the gravitational model is 
selected with the sole purpose of supporting the metrics of interest and give rise to the de- 
sired properties. For instance, the Lifshitz metric ( II. ip is a solution of AdS gravity coupled to 
a massive vector [7] if one tunes the mass and the cosmological constant to m 2 = z (d — 2) L~ 2 
and 2A = — [z 2 + (d — 3)z + (d — 2) 2 ] L~~ 2 respectively. Notice that both the dynamical ex- 
ponent and the dimensionality of the spacetime enter in these expressions. 

The flexibility of the phenomenological models is certainly an advantage and much in- 
sight can be gained with this modus operandi. Nonetheless, gauge/gravity dualities are 
firmly established just in the context of string theories and one could be skeptical about the 
legitimacy of bottom-up constructions. In other words, it is important to verify that the 
relevant solutions can also be obtained within string or M-theory. Unfortunately, as one can 
presume, this is a more intricate task. 

However, many important steps have been taken in this direction. By now we have at 
our disposal embeddings into string theory of holographic superconductors [HI [9] as well as 
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the Schrodinger spaces mentioned above [TOl [TTl [T2] . Lifshitz geometries have proven to be 
more elusive, but they were recently found in [131 HS1 HI HZ] • 

Consistent truncations have been instrumental in obtaining various of these solutions. 
They provide lower- dimensional simplified setups with respect to 10 or 11-dimensional su- 
pergravities. Compared to phenomenological models they have the virtue that the embedding 
of the solution into string theory is automatic. The drawback is rigidity, in the sense that the 
parameters of the theory, such as masses and charges, are fixed during the reduction. Com- 
ing back to the example of Lifshitz this means that, either one is working with a consistent 
truncation containing a massive vector and a cosmological constant with the exact values 
given above, or the solution is not there. One could then rephrase the search for Lifshitz 
solutions in string theory as the chase for a vector with the correct mass. 

An important point was the observation that one can keep the appropriate massive modes 
in the reduction of higher- dimensional supergravities [10J. Stimulated by this result, several 
consistent truncations containing massive vectors were subsequently found [TS1 ESI EDI [2D 
1221 [221 EES] • Nevertheless, none of the vectors included in those truncations seemed to have 
the correct mass (except one in [16] that gives z ~ 39). A proper vector to support Lifshitz 
metrics with z = 2 was identified in [13] and further exploited to construct infinite families 
of solutions in [TH [T7] . 

On the other hand, it is known since the dawn of the AdS/CFT duality that introducing 
temperature in the gauge theory corresponds to placing a black hole in the bulk of the gravity 
side. When the gauge theory is conformal, the corresponding black hole has AdS asymptotics 
in the boundary. Similarly, since the Lifshitz field theory is endowed with anisotropic scale 
invariance, the pertinent black holes describing its finite temperature generalization must 
have the metric (11. ip as asymptotic geometry. 

Black holes possessing this property have been extensively studied in phenomenological 
models, starting with the d = 4, z = 2 numerical ones of [21]. Extensions to other dimen- 
sionalities and critical exponents are available as well [23 [261 [23 121] ■ The charged case 
is also interesting for condensed matter applications [211 EDI EI]- Lifshitz black holes are 
also solutions of higher-order theories of gravity [221 E31 EH E2] and Brans-Dicke models 
[SB] . Even if most of the results are numeric, a few exact black holes are known [37, 38J. 
Thermodynamic properties of these systems are detailed in [39, 40J. 

Despite all this considerable effort, a proper embedding of Lifshitz black holes into string 
theory was still missing^]. In this work, we construct an infinite family of d = 4, z = 2 
Lifshitz black holes in type IIB supergravity. Building on the results of [HI [17] , we present a 
consistent truncation based on E5 x S 1 , where E5 is an arbitrary 5-dimensional Einstein space, 
that contains a massive vector suitable for supporting the solution. In addition it incorporates 
two scalars, the dilaton and the modulus of the circle, whose presence complicates the study 

1 The solutions of [41] are rather naked singularities than genuine black holes. 
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of the equations but triggers interesting new features with respect to the more minimal 
phenomenological models. In particular, we find no extremal limit and an instability for 
small black holes similar to the usual AdS-Schwarzschild one. This signals to a possible 
Hawking-Page like transition |42] in these systems. 

The structure of the paper is simple: in section [2] we construct the consistent truncation 
and present the 4-dimensional model and its equations of motion, in section |3] we solve 
numerically the equations and comment on the properties of the black holes found and in 
section 0] we conclude. 



In accordance with the results of |17j . a 4-dimensional Lifshitz space with z = 2 can be 
attained starting from a model in five dimensions containing an axion field and admitting an 
AdSs vacuum. The mechanism is the following: if the 5d theory is reduced on S 1 while using 
the axion field strength to thread a flux, the vector gauging the circle isometry has the correct 
mass to support the 4d Lifshitz solution. In this section we will provide a very simple model, 
coming from a consistent truncation of type IIB supergravity, that meets these requirements 
and thus is suitable for the embedding of Lifshitz-like solutions into string theory. Using an 
ansatz that generalizes the metric (11. ip to allow for a horizon, we derive the complete set of 
equations to be solved. 

A straightforward way to get an AdSs solution in type IIB is to compactify on a 5- 
dimensional Einstein space E 5 with a Freund-Rubin ansatz for the RR 5-form. On the other 
hand, the most natural candidate for the axion would be the RR scalar Co- Indeed, as 
already mentioned in [17], the model obtained in this manner is a consistent truncation of 
type IIB supergravity as long as one retains in addition the dilaton 0. In detail, we take the 
10-dimensional metric to be 



where M5 is for the time being an arbitrary 5-dimensional spacetime and the internal Einstein 
metric is normalized as Rij = ^g%j. By assumption the axio-dilaton has sole dependence on 
the M5 coordinates and from the set of forms the only non- vanishing one is 



ensuring self duality. We denoted Vol(Es) the volume form on the Einstein space. Within this 
ansatz the Bianchi identities are satisfied (taking F\ = dCo), as well as all the 10-dimensional 
equations of motion except the one for the axio-dilaton and the external components of 
Einstein's equations, that follow from the 5d action 



2 The model 



ds 2 = ds 2 (M 5 ) + ds\E 5 ) 



(2.1) 



F 5 = 4 (1 + * 10 ) Vol(E 5 ) , 



(2.2) 



S 




1 



/ 



R * 1 d(j) A*d<p e 2lf> dC A *dC + 12 * 1 



(2.3) 
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Notice the presence of a negative cosmological constant A = —6 that emerges from a com- 
bination of the internal space curvature and the 5-form flux. It is then apparent that the 
model contains an AdSs vacuum for constant scalars. In order to find the 4-dimensional 
Lifshitz solution one needs to further reduce the theory on a circle with metric 

ds 2 (M 5 ) = e~ T ds 2 (M 4 ) + e 2T (d& + A) ® (d& + A) . (2.4) 

The scalar T(x) is a modulus parameterizing the size of the S 1 while the vector A(x) is 
gauging the reparameterization invariance of the coordinate d on the circle. Both depend 
exclusively on the coordinates of the 4d spacetime, denoted collectively x. This vector turns 
out to support a Lifshitz metric if one Stiickelberg couples it to the axion via a flux on 
the S 1 . In doing so, the correct ansatz for the axion is Co(x,$) = Cq{x) + n$ with n a 
constant. We should point out that, even if Co has a explicit linear dependence on the circle 
coordinate, its axionic nature ensures that it necessarily enters into the action covered by a 
derivative and thus the circle reduction is consistent. Moreover, in this way one generates a 
covariant derivative for it and consequently a mass term for the vector 



dC Q {x,d) = dC (x) + ndd= DC {x) + n (d& + A) 



(2.5) 



with DCo = dCo—n A. The second term will enter as a modification to the potential. Assum- 
ing that the dilaton depends just on the M4 coordinates the reduction is then straightforward 
and the resulting action reads 



S 



2k\ 



3 1 

R*l — -dT A *dT d(f) A *c 

2 2 



- -e 2<t> DC A *DC 



-- e 3T dA A *dA + ( 12 e- T - — e 2 *~ 3T 
2 V 2 



* 1 



(2.6) 



where the gravitational constant has been redefined to K5 = k\ J gl d$. For our purposes it 
is more convenient to dualize the axion to a 2-form and work with a system more similar to 
the one in [6] . Let us add to the action the following piece 

1 



S' 



2k 2 



(Gi + nA) A dA 2 , 



(2.7) 



where we denoted G\ = DCq. Variation with respect to A 2 gives the Bianchi dG\ = —nT as 
expected, while varying w.r.t. G\ one deduces instead the duality condition G\ = e~ 2 ^ *dA 2 
that once substituted into the action yields 



S 



1 

2^1 



R*l 



-dT A *dT 



1 

2 ( 



A *d(p — 



- 2 *dA 2 A*dA 2 



1 



e 3T FA*J c + I 12 e 



(l2 e- T - y e 2 ^ 3T ^) * 1 + n A 2 A T 



{21 
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Leaving apart the scalars <p and T, this is the form of the action used in [6]. These two 
modes are spurious for the purposes of obtaining Lifshitz solutions, but are unavoidable for 
consistency of the truncation. Indeed, the scalar T will be present in any model that uses 
A as a support for the Lifshitz metric, since it is well known that the modulus of a circle 
cannot be consistently switched off while retaining the vector gauging its isometry. Similarly, 
we cannot truncate away the dilaton and simultaneously keep the axion. The presence of 
these spectator fields complicates the analysis of the equations of motion with respect to the 
phenomenological models of [HI 121] • On the other hand, some new appealing features of our 
black holes will be caused by these modes, as we will detail. 



We want to genera 
with z = 2, so we piclj 2 



ize the Lifshitz metric (11. ip above to allow for black hole solutions 



ds 2 = -r 4 f(r) dt 2 + r 2 dQ 2 2 + ^ dr 2 , (2.9) 

where the transverse piece reads 

{dip 2 + sin 2 ip d\ 2 if k = 1 

dip 2 + ip 2 d X 2 if k = (2.10) 

dip 2 + sinh 2 ip dx 2 if k = -1 

and k is the curvature corresponding to spherical, flat or hyperbolic horizon, respectively. 
Since we are interested in geometries that are asymptotically Lifshitz we will impose the 
boundary conditions f(r), g(r) — > 1 as r — > oo. 

To support this metric we take the following ansatz for the forms (with the obvious choice 
of vielbeins) 

T = 2 h(r) e r A e l , 

F 3 = - j(r) 9 r A9^ A 9 X (2.11) 

n 

and permit a radial dependence of the scalars 



n 2 



T = T(r) and e" 2 * = — e~ 2lp{r) . (2.12) 

Substituting this ansatz into the matter e.o.m. we get the following system of coupled 



2 Here we set the size of the space to L = 1 for simplicity. All the dimensionful quantities in the following 
are measured in units of L. 
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differential equations 



[e 3T hr 2 ]' - 2jgr = 0, 
[e-^jfr 2 ]'-2hfgr = 



ip' r f 
9 . 

T'r 5 f 
9 



fgr 3 

1 

f gr 3 



3T 



+ 2 e 3T /i 2 - 4 e" T + 2 e 2</, - 3T = 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



where ' indicates derivative with respect to the radial coordinate. There are two combinations 
of Ricci tensor components, namely R tt + R rr ± 2-R^, that give first order equations once 
substituted into Einstein's equations. They read 



^ = -5 + | (T' r) 2 + 1 (y/ r) 2 + / (V 2 ^ j 2 - e 3T /> 2 + 6 e~ T - e 2 ^ T + 



(2.17) 



and 

2rg' 
9 



3 + | ( T ' r) 2 + - A {ip 1 rf + ^ 2 (e~ 2 * j 2 + e 3T /i 2 - 6 e~ T + e 2 ^ 3T - ^ . (2. 18) 



The third component is not independent and follows from these two. 

One important point is that, since /' oc /, we can use (I2.17P in the e.o.m. of the matter 
fields to decouple / from the system. We can then solve for T, ip, g, h and j and substitute 
in the equation for /. We have not found exact black hole solutions of the complete system 

It is simple to check that the Lifshitz metric (11 .ip with z = 2 and d = 4 verifies the 
equations above. This amounts to take 



/ = 9 



J =1 



and 



T = 0. 



(2.19) 



Therefore, the solutions we are interested in will approach these values asymptotically in the 
boundary located at r — > oo. 



3 The Lifshitz black hole solutions 

In order to find asymptotically Lifshitz black hole solutions we will resort to numerics. But 
before integrating the system, let us study the asymptotic behaviour of such solutions. 

3 Plugging the metric and form functions of the topological black hole in [37] one can not verify the 
equation of the dilaton (12.15[) . 
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3.1 Asymptotic behaviour 

In the UV regime, r — > oo, we are looking for solutions that asymptote to Lifshitz geometry 
with scaling z — 2, i.e. /, g, h, j — > 1 and if, T — > near the boundary. Therefore we can 
linearize the equations of motion around the fixed point and study the subleading behaviour 
of the fields in the asymptotically large r region. For convenience let us define new scalar 
fields 



T r 



j(r) 



(3.1) 



that approach t, 7 — > 1 for the Lifshitz solution. If we expand the fields around the scale 
invariant solution as ip(r) = ipuf + Sip(r), ft * s straightforward to solve analytically the 
resulting equations of motion up to linear order in Sip. The most important feature of the 
linearized system is the presence of a zero mode solution, i.e. independent of the radial 
coordinate. Of course, coming back to the non-linearized problem, such zero mode develops 
an r-dependent profile. If one gets close enough to the boundary, the zero mode gives the 
dominant subleading contribution to the fields, since one expects the other modes to have 
already decayed. However, even if the zero mode has negative amplitude, its decay will be 
very slow and so will be the convergence to Lifshitz geometrjfl. For that reason we will seek 
solutions for which the amplitude of the zero mode vanishes and the decaying modes of the 
fields behave as 1/r 2 and log r/r 2 . In particular, the leading orders of the solutions to the 
non-linearized system in the UV, imposing Lifshitz asymptotics and a vanishing amplitude 
for the zero mode, are given by 



f(r 
g(r 
h(r 

j(r 
7 ( r 

r(r 



33 + a — k 8a — k 



1 - 
1 - 

1 - 



6/3 



2r 2 
- 2a 



a: 
33 



+ . . . 

— 5a 



+ 



4r 2 

8a — k 



4r 2 



log r + . . . , 
log r + . . . 



(3.2) 



2r 2 



la — k 



log r + . . . 



a. 



sa 



6r 2 



logr + 



where a and 3 are arbitrary constants and the dots represent terms that decay faster than 
1/r 2 . 

In the deep IR we want the above solutions to flow to a black hole like geometry. A 
non-extremal black hole is characterized by the presence of a non-degenerate horizon at a 
finite r = rn- In our coordinates, this implies that the g u component of the metric has a 



4 It is easy to check that our non-linearized system is solved by a series expansion in a marginally decaying 
mode of the form c(r) = log (log r) / log r . 
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simple zero at this point, whereas the g rr component has a simple pole there. Regularity of 
the solution further demands that the 2- form field strength j(r) vanishes at the horizon as 
~ l/g(r). We will assume that the vector field strength and the scalar fields have a finite 
value at the horizon. Under these considerations, the near-horizon expansion of the fields 
takes the form 



g(r 

j(r 
h(r 
7 ( r 

r(r 



Vr - r H (fo +fi(r- r H ) + / a (r - r H ) 2 + • • •) 
(go + 9i(r - r H ) + g 2 (r - r H ) 2 + . . . ) 



Vr - r H (jo + ji(r - r H ) + h{r - r H ) 2 + . . . ) 
/i + h(r - r H ) + /i 2 (»" - ?h) 2 + • • • , 
7o + 7i( r - r H) + T2(^ - r H) 2 + • • • , 
T + n(r - r H ) + r 2 (r - r H ) 2 + . . . . 



(3.3) 



Plugging that expansion into the equations of motion and solving order by order we find that 
/o just appears as a global factor in the f(r) field. This is because its equation of motion 
is linear in f(r) so the normalization of the field is not fixed. This fact will be important 
for the numerical integration. Once we rescale f(r) by /o, all the coefficients can be written 
only in terms of ho, 70 and r , being the first ones in the series given by 



9o 



Jo 



k j 6_ 7n L2_3 

^ + ro ^ a T 



2/lQ70 



k _i 6_ 7o L2_3 



(3.4) 
(3.5) 



It is clear that the horizon values of the vector field strength and the scalars are constrained 
by demanding the solutions to be real. For the planar and spherically symmetric black hole 
solution^! it implies the following upper bounds 

l^ol < —o\ — 5^ + 6 — and — < \ + 6 , (3.6) 

To V ^ T r « V r H 

where 70 and r are both positive and from now on we assume h to be so. 
5 We will not further consider the hyperbolic case k = — 1. 
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3.2 Numerical integration 

Now we can proceed to the numerical integration of the system. We use the near-horizon 
expansions in (13.31) to generate initial data close to the horizon and then integrate the 
equations of motion in (I2.13p -f l2.18p towards the UV demanding that the solution has Lifshitz 
behaviour at the boundary. Since we have a large number of free parameters the strategy 
is the following. As we already know, /o is nothing but a normalization scale for f(r). 
Therefore we can set /o = 1 and afterwards rescale with the asymptotic boundary value 
of the numerical solution, automatically fixing the normalization of / in such a way that 
/(r) — >■ 1 as r — >• oo. We are then left with four parameters, i. e. the horizon radius, the 
vector field strength and the scalar horizon values. In each integration we keep fixed the 
black hole size rn and the modulus To, while tuning the dilaton 70 and the 1-form flux ho in 
order to get Lifshitz solutions. The precision of our numerical results is — 1| < 10~ 4 at the 
cutoff r = 10 6 , where ip represents all the fields but /. We then explore the parameter space 
by repeating the procedure for different pairs of initial data {?"h, To}. In principle we could 
have picked any other free parameter rather than r to generate the initial data. However 
this is a natural choice since the modulus T will be present in any Lifshitz solution based on 
the massive vector A [HI [T7] . 

In figures [U [2] and [3] are shown the metric, form and scalar radial profiles of the solutions 
for two different black hole sizes and each for two different values of the modulus. It is clear 
that convergence to Lifshitz geometry is better for large values of the scalar tq at the horizon. 
Convergence is also better the larger the black hole for a given To. It turns out that 7(7™) < 1 
for any pair of initial values of the parameters. This implies that the dilaton in the black 
hole solution is always smaller than in the pure Lifshitz space, i. e. </>bh(?") < <fiuf- On the 
other hand, the constraint found on the modulus is simply given by r > 70 > (see figure 
SJ), so the original modulus T(r) can take both positive and negative values in the black hole 
solution. The peak and dip developed by the metric functions of the black holes close to the 
horizon grow when the modulus To decreases. 

From this behaviour and the fact that convergence is worse the smaller is To, one could 
expect that for a given black hole size, there is a minimum value of that scalar field at 
the horizon for which the field solutions are able to stabilize to their Lifshitz limit when 
flowing to the UV. In fact this is the case. In figure 3] it is plotted the 1-form flux as a 
function of r )o/ T o f° r fixed values of the horizon radiusj. Very large values of r correspond 
to the beginning of the curves, where ho — > 0, whereas the ending of the curves corresponds 



6 In principle, for a given horizon radius, kg is a function of both 70 and To and not just of its ratio. 
However, in the k = case the system is invariant under rescaling of the radial coordinate and solutions 
are independent of th- Furthermore, solutions are related by rescaling of the fields in such a way that the 
space of free parameters is reduced to {ho, 7o/t }. On the other hand, the flat case k — is equivalent to 
the rn — > 00 limit of the spherically symmetric black hole. Comparison of the flat and spherical black hole 
solutions is possible by plotting ho vs. 7o/tq. 
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Figure 1: Metric functions f(r) (solid) and g{r) (dashdotted) for various different pairs of 
initial values. The red (inner) and the blue (outer) lines correspond to r = 2 and t$ = 1 for 
r H = 1/2 (left), and to r = 2 and r = 1/2 for r H = 20 (right), respectively. 
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o r J ' 1 

5 10 15 20 5 10 15 20 

r/r H r/r H 

Figure 2: Form fields h(r) (solid) and j(r) (dashdotted) for various different pairs of initial 
values. The red (inner) and the blue (outer) lines correspond to tq = 2 and tq = 1 for 
rn = 1/2 (left), and to r = 2 and r = 1/2 for rn = 20 (right), respectively. 



1.5 - 




1 1 1 

5 10 15 20 5 10 15 20 

r/r H r/r H 



Figure 3: Scalar fields 7(7") (solid) and r(r) (dashdotted) for various different pairs of initial 
values. The red (upper) and the blue (lower) lines correspond to r = 2 and r = 1 for 
r H = 1/2 (left), and to r = 2 and r = 1/2 for r H = 20 (right), respectively. 
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to the minimum value of Tq for each rn- At this minimum, 70 is also a minimum, while ho 
reaches its maximum value. For the flat black hole solution, k — 0, the limiting configuration 
corresponds to the maximum value of the ratio of the scalars, 7o/to = 0.8164, and to the 
minimum value that the scalars can reach. 

The bounds (13. 6p are satisfied along the whole curves. Moreover, the bound for ho is 
saturated only in the tq — > 00 limit, where ho — > 0. For very large values of the modulus tq, 
the 7 field is roughly constant, 7(7") < 1, so there is no running of the dilaton with the RG 
flow, and the horizon 1-form flux almost vanishes. Therefore the black hole is only supported 
by the modulus T. Notice that this behaviour is independent of the black hole size and that 
for non- vanishing rn extremality is only achieved in the limit To — > 00, that of course is not 
a regular solution at the horizon. For regular solutions extremal black holes could appear in 
the limit of vanishing black hole size. This kind of zero temperature solutions were found in 
phenomenological models j2U [37] . However, it will not be the case in the stringy model at 
hand due to the presence of the additional scalars. 

100 

80 
60 

o 

40 



20 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

7(1 /tq 

Figure 4: Horizon 1-form flux h as a function of the ratio of the scalar fields 7o/to for fixed 
rn = 1/2, 1, 3/2, 2, 5, 20 and rn = oo, from left to right. 

For a given initial To there is a minimum black hole size that can support the configuration 
of fluxes that make the solution flow to Lifshitz asymptotics. In figure [5] we show the horizon 
values 70 and h as a function of r H for fixed values of r . For large r H they take the 
asymptotic values given by the flat solution, k = 0, and both of them remain barely constant 
in a wide range of t"h- When the horizon radius gets close to its minimum size, the dilaton 70 
rapidly decreases to its minimum value, while the flux ho either grows or decreases depending 
on the value of r . There is a second branch of solutions that is not connected to the flat 
solution and that merges smoothly with the one described above at the minimum r^. Over 
this second branch, 70 decreases and ho increases while the black hole grows. In principle, 
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Figure 5: (Left) Horizon dilaton 7 as a function of rn for fixed r = 6, 2, 1, 1/2, from top 
to down. (Right) Horizon 1-form flux ho as a function of rn for fixed tq = 6, 2, 1, 1/2, 
from down to top, normalized by the flat solution value: = 0.0549, 0.459, 1.57, 4.79, 
respectively. First branch (solid) is connected to the k = case, while the second branch 
(dashdotted) is unstable. 



we expect 7 to asymptotically vanish for large black hole radius. Convergence to Lifshitz 
geometry of these solutions is very poor when one moves further from the minimum size, 
contrasting with the first branch, where convergence is better the larger the black hole. The 
second branch will be shown to be thermodynamically unstable. 

3.3 Thermodynamics of Lifshitz black holes 

From the ansatz ( 12. 91) and using the near-horizon expansion ( 13.3f) . we can determine various 
thermodynamical quantities that characterize the black holes. The Hawking temperature 
and entropy density are given by 

T H = (3.7) 

471^0 

s = nr 2 H , (3.8) 

both for k = 0,1, after appropriate volume normalization and where we set the Newton 
constant to one. The conditions ( 13. 6 p have to be satisfied in order to have a positive real 
temperature. In fact, if the bound for ho is saturated the black hole becomes extremal. Due 
to the intrinsic dependence of the horizon values of the fields on the black hole size it is 
not possible to write the explicit dependence of the entropy density on the temperature. 
However approximate relations can be found in certain regions of the parameter space. 

The thermodynamic behaviour of the Lifshitz black holes above can be easily studied. 
Figure shows the Hawking temperature as a function of the black hole size and the entropy 
density as a function of the temperature. There are significant differences with previous 
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Figure 6: (Left) Temperature as a function of rn and (right) entropy density as a function of 
temperature of the first (solid) and second (dashdotted) branches for fixed r = 6, 2, 1, 1/2, 
from down to top (left) and from left to right (right). It is also shown the asymptotic 
behaviour for large values of radius (dotted). 

phenomenological models. As already mentioned, in [2U [37] extremal black holes were 
found in the limit of vanishing horizon radius since temperature for such models turned 
out to be a monotonic function of the black hole size. In our string theory model, the 
presence of additional matter fields prevents the horizon to shrink to rn = and consequently 
extremality is never reached. Temperature for large black holes increases with the size, while 
for small black holes it develops a negative slope. Having negative specific heat, these small 
black holes are always thermodynamically disfavoured. Notice that the concept of smallness 
is T -dependent. This is precisely the behaviour shown by AdS-Schwarzschild black holes. 
In analogy with this well known case [12] > we expect to have a Hawking-Page like phase 
transition between large Lifshitz black holes at high temperature and 'thermal' Lifshitz 
(pure Lifshitz space with compact Euclidean time) at low temperature. 

For very large black holes the temperature is approximately T H cx r H . Since the Bekenstein- 
Hawking entropy density is simply given by the area of the horizon, it is possible to write 
down its temperature dependence in the rn 3> 1 limit as 

s = a T0 Th , (3.9) 

where the proportionality constant can be computed numerically. For the solutions presented 
in figure El it reads 

a 6 = 345.2, a 2 = 36.53, an = 8.055 , a 1/2 = 1.619. (3.10) 

Hence, large Lifshitz black holes satisfy the expected scaling behaviour for a field theory 
with dynamical exponent z = 2 in flat 2 + 1-dimensions, s oc T 2 / 2 . 

In order to ensure that the small black holes and the second branch of solutions are 
unstable compared to the large ones we can compute the free energy of the system. From 
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Figure 7: Free energy density as a function of temperature of the first (solid) and second 
(dashdotted) branches for fixed r = 6, 2, 1, 1/2, from left to right. 

the first law of thermodynamics, e + P = sT, and the relation between energy density and 
pressure for the dual field theory, e = P, the free energy is easily computed as 



This thermodynamic relation has been proved to hold by computing the free energy from 
the Euclidean bulk action in [281 E3 EQJ H3] for various phenomenological model^l Figure 
[7] shows the free energy as a function of the temperature for various initial values of r . 
For temperatures above the minimum temperature, there exist two different solutions. The 
lower one, hence the stable, corresponds to large black holes, whereas the upper one is clearly 
unstable and corresponds both to small black holes and to the second branch of solutions, as 
expected from previous considerations. For temperatures below the minimum temperature, 
Lifshitz black holes can not exist and presumably a phase transition to 'thermal' Lifshitz 
will take place. 

4 Conclusions 

In this work we have constructed infinite families of asymptotically Lifshitz black holes 
in string theory. We made use of a consistent truncation of type IIB supergravity on an 
arbitrary 5-dimensional Einstein space that keeps the axio-dilaton. By further reducing on a 
circle we obtain an additional scalar (the modulus of the S 1 ) and a vector that is appropriate 
to support Lifshitz black holes with dynamical exponent z = 2. The model obtained in this 

7 In our case we have extra scalar fields not considered in the mentioned works. However it is not expected 
that the presence of additional scalars spoils this result once the appropriate counterterms are taken into 
account. 



Q = -P = —sT. 



(3.11) 
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way is not minimal, in the sense that previous phenomenological studies were lacking the 
dilaton and the modulus that are mandatory for the consistency of our truncation. 

These extra modes contribute in a very non-trivial way to the equations and yield new 
interesting features. We found no extremal black holes. In phenomenological models [2U [37] 
the extremal limit corresponds to vanishing black hole horizon, but here we observe a lower 
bound on the size. On top of that, small black holes turn out to be unstable. For large black 
holes, temperature increases with the horizon radius, while for small black holes it develops 
a negative slope, pointing to an instability that we confirmed by computing the free energy. 
The fact that black holes have a minimum temperature leads us to expect a Hawking-Page 
like phase transition, resembling the AdS-Schwarzschild case. 

It would be interesting to know if this behaviour is general in string theory constructions. 
It may happen that this property is common to solutions based on the vector we kept, that 
comes always accompanied by the modulus of the circle. This could be checked by using 
truncations with matter content different from the axio-dilaton, for instance the one based 
on the T 1 ' 1 of [2]. Other suitable truncations should be easy to find [T7] . 

A pressing question is the nature of the expected phase transition. The first step would 
be to compute the free energy of the thermal Lifshitz space, in order to proof its existence. 
We must point out that for Lifshitz geometries the gauge/gravity holographic dictionary 
is underdeveloped, though some progress has been made |44j . This makes involved the 
identification of the relevant degrees of freedom in each of the phases. In consequence, 
it might be premature to interpret the phase transition as the analogue of confinement- 
deconfinement [15] in non-relativistic systems. We hope to come back to these issues in the 
near future. 

It is certainly possible to extend our results to charged Lifshitz black holes. When the 
internal manifold is not only Einstein but admits a Sasakian structure, there is a consistent 
truncation that includes an additional vector. Switching on this vector introduces a charge 
density and thus charged black hole solutions similar to Reissner-Nordstrom can be obtained. 
Interestingly, the truncation can be enhanced to include a charged scalar [HI [XHl LTT] . Lifshitz 
black hole solutions where this scalar condenses and breaks the Abelian symmetry would be 
dual to the superfluid phase of a non-relativistic superconductor [29|. 
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